Diffuse transport and spin accumulation in a Rashba two-dimensional electron gas 
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The Rashba Hamiltonian describes the splitting of the conduction band as a result of spin-orbit 
coupling in the presence of an asymmetric confinement potential and is commonly used to model 
the electronic structure of confined narrow-gap semiconductors. Due to the mixing of spin states 
some care has to be exercised in the calculation of transport properties. We derive the diffusive 
conductance tensor for a disordered two-dimensional electron gas with spin-orbit interaction and 
show that the applied bias induces a spin accumulation, but that the electric current is not spin- 
I , polarized. 

PACS numbers: 72.10.-d,72.20.-i,72.25.-b 

Utilizing the spin degree of freedom for electronic applications is a declared goal of the research field of mag- 
netoelectronics or r&Bin-electronicsEl. Devices made from metallic layered systems displaying the giantcl and tunnel 
magnetoresistanceElcl have been proven useful for read-head sensors and magnetic random access memories. Inte- 
gration of such devices with semiconductor electronics is desirable but turned out to be difficult because a large 
resistivity mismatch between magnetic and normal materials is detrimental to spin inject ionu. Still, th^pxablem 
can been solved in various ways and spin injection into bulk semiconductors has indeed been reportecloLnau'EJEil. 
Electrical spin injectiou_into a high-mobility two-dimensional electron gas (2DEG) and its detection appears to be 
much more demandingt3. In this context it would be attractive if application of an electric field alone would suffice 
to induce a non-equilibrium magnetization or spin accumulation in the presence of the spin-orbit interaction. Such 
an effect, dubbed the "kinetic magnetoelectric effect" was actually predicted in seminal theoretical work by Levitov 
et al.t3, and is caused by the combined action of the spin-orbit interaction, absence of inversion symmetry, and the 
time-reversal symmetry-breaking by an electric field in disordered systems. In asymmetric heterostmctures made from 
narrow-gap semiconductors the spin-orbit-ipteraction is dominated by the rS<p-called Rashba termcj, which has a very 
simpk-jStructure, can be quite significanffi-3, and modulated by gate fieldai3. Recent observations of a spin-galvanic 
effectEZI and spin-orbit scattering induced localization/antilocalization transition in 2DEGsE§E3 reflect the interest and 
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importance_af the topic. 

EdelsteinEl showed that an applied field induces an in-plane magnetization in a Rashba 2DEG. Although he did 
not make any suggestions in this direction, the interpretation of the spin accumulation being caused by an effective 
magnetic field has lead subsequently to the misconception that the current is also spin polarized. A microscopic 
calculation of both spin accumulation and current (or conductivity tensor) on an equal footing is thus required. 
Furthermore, the existence of the spin-^r bit.. induced spin-accumulation and the conditions for its observability have 
recently been a matter of controversyuO'EaESEj. 
^ ' In this Report we carry out microscopic model calculations of the conductivity tensor and spin accumulation for 
J> . a disordered Rashba 2DEG in the linear response regime. This task is complicated by the correction to the electric 
field vertex, which does not vanish even for short-range isotropic scatterers. We confirm that a spin-accumulation 
normal to the applied electric field vector is excited. However, the electric current is not spin polarized, thus solving 
the controversies mentioned above. We furthermore show that the mobility increases quadratically with the Rashba 
spin-orbit interaction. 

The Rashba Hamiltonian in the momentum representation and Pauli spin space reads 

( £~k 2 i(aE z )k_\ M . 
V -t(aE sc )k+ J 

where k± = k x ± ik y with k = (fc x , k y ) the electron momentum in the 2DEG plane. (aE z ) parametrizes the spin-orbit 
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coupling and is experimentally accessibleB. The eigenfunctions of the Hamiltonian are 
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e ik TV )- (2) 



with s — ± , k — ^Jk 2 + fc 2 anc l L 2 the area of the 2DEG and corresponding eigenvalues are given as Ek s = fi 2 k 2 /2m+s 
(aE z )k. The current operator in Pauli spin space is given as 



IX — CUi 



J, = e^ = e( 6 ^ 6 fcJ, (4) 

with b — h/m and A = (aE z )/%. In the space of the eigenfunctions of Hq, referred to hereafter as the s-space (s = ±), 
the current operators are transformed as 

My) = U^j x(y) U = e (bk x(y) l + \-^a z - {+)X^a}j , (5) 

by the unitary matrix 

i / .fe_ .fc_ \ 

We also need the transformed spin matrices, U'a x U, U*a y U and U'cr z U to evaluate the spin accumulation. 

The standard model for disorder consists of randomly distributed, identical point defects, which are neither spin- 
dependent nor flip the spin: 

V(r) = V1^6(r--Ri). (7) 

i 

In the following, we expand the Green function G = (zl — H) -1 , where H — Hq + V(r) and z = e± ir/, in terms of the 
unperturbed Green function, Go = [zl — Hq] , with matrix elements gk± = l/(z — Ek±) in s-space, and calculate 
the self-energy S in the Born approximation. 

After ensemble averaging over the impurity distribution, denoted by (• ■ -)aVi an d disregarding a trivial constant 
term, the self-energy in the Born approximation reads: 

nV 2 / k , V V k \ 

«k S |vG y|k'v')w = ^j;^ ( 1 + ss " + ss '^kr + s ' s "^r) ^ 

k's' ^ ' 

Swdss" y^fffc's' = S<$kk"<W', (9) 
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where n = N/L 2 is the density of impurities per unit area. Eq. (9) follows from the odd symmetry of k + k'_ and 
k' + k_ with respect to k' x or fc^. The Green function is therefore given as 

((ks\G\k" s ")) AV = _ 1 1 £ kk »<W = G fcs . (10) 

We find both G and S to be diagonal in k and s. By direct inspection it can be seen that (VGqVGqVGqV) is 
diagonal, meaning that the exact self-energy must be diagonal as well. 
The (longitudinal) conductivity is given by the Kubo formula as 

^Tr (J X G A J X G R + J x G R J x G A ) Ay , (11) 
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where the superscripts A and R denote Advanced and Retarded, respectively, and will be omitted below for brevity. 
We evaluate (J x GJ x G)av — Jx(GJ x G) av = JxK in the ladder approximation, which obeys the Ward relation with 
the self-energy in the Born approximation: 

K ~ GJ X G + G(VKV)avG. (12) 

The matrix elements of (VKV)av ar e 



((ks\VKV\k's')) A v = EE (aZ*) E(^ kl) - Rl ^ (k " k2) ^ 



AV 



SSi ^r) ( x + S2S '^yr) <^^]^i^ S2 >. us) 

To evaluate the expression for (V-KV^v, we first use G{k\S\){kiSx\ J x \k2S2)G{k2S2) for (kiS\\K\k2S2)- Because J x is 
diagonal in k, ki = k2, and the average of the exponential factor leads to k = k'. Repeating this procedure iteratively, 
we find that, like G, K is diagonal in k. The matrix elements (ks\V KV\ks') may be evaluated iteratively. We call 
(ks\K \ks')^> = (ks\GJ x G\ks') and note that (ks\GJ x G\ks') cx k x when s — s' and cx k y when s ^ s' and keeping 
terms which are even functions of k x and k y in the summation of the equation. Then, by direct inspection, 

(ks\VKV\ks') AV = ^- (k x a z - k y a y ) , (14) 

where 

e\' = ^Y.\{ k ^ k + \G J *G\k+)~ik y {k + \GJ x G\k-) 
k 

+ ik y (k - \GJ x G\k+) -k x (k- \GJ x G\k-)^ . (15) 

Comparing Eq. (14) with Eq. jl^ ) and the definition of J x , Eq. (4), we find that K = GJ X G has the same structure 
as GJ X G, where A in J x is replaced by A + A' in J x . By replacing GJ X G in Eq. (|l5| ) with K = GJ X G, we obtain a 
closed equation for A': 

A ' = i^E{^ + ( A + A 'w- ( 16 ) 

fe 

with S\ = sGksGks and So = ^2 ss i GksGks' ■ Here, we have used relations k x = k y = k 2 /2 in the summation over 
k. 

The conductivity now reads 

TtJ x GJ x G, (17) 
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where the tilde indicates substitution of A by A + A', and J X GJ X G is a 2x2 matrix expressed in s-space. By carrying 
out trace, the conductivity follows as a xx = a xx + a xx with 



xx AttL 2 

k 



^ = ^E{( &± ^) (fe±^^)fc 2 G fc± G fe± + A(A + A')G fcT G fe± |. (18) 



By using the unitary matrix U, the matrix representation of the conductivity in s-space can be transformed into that 
in the original Pauli spin-space as a xx = hU J X GJ X GW /2ttL 2 , with a'^x'^'xi e ^ c - ^v taking the spin trace of a xx , 
the relation a xx + a xx = &l x + a**, follows naturally. We observe that a xx = and that the non-diagonal elements 
of the conductance tensor in the original spin-space a\\ vanish identically by parity. These results prove that the 
current excited by the electric field is not spin polarized. The expression for a yy can be derived analogously. We also 
find that the non-diagonal (Hall) conductivity a xy vanishes by symmetry. 

In calculating the vertex correction and conductivity we encounter integrals over the momentum, which may easily 
evaluated by an approximation in which the life time broadening of the density of states is neglected, 

G% s Gl = ^5(e-E ks ), (19) 
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where the life time r is defined £ = — isign(n)?i/2T or r = Ti/2-KnV 2 D with the 2DEG density of states per spin 
D = m/2-Kh 2 . We then obtain A' = — A, and 

2e 2 nnT o „ 

o** = — + 2e 2 rDX 2 , 20 

m 

where no is the number of electrons per spin. The conventional Drude conductivity (first term) is increased by the 
spin-orbit interaction. Since the sign of the coupling constant A is irrelevant, the enhancement term must be of even 
order in A . The result A' = —A shows that the vertex function J x is diagonal in spin spacec3. In the case of A = 0, 
the vertex correction in the ladder approximation vanishes identically due to the isotropic scattering. 

It is important to distinguish the spin-polarized currents computed above from the spin-accumulation which is 
excited by the applied field E, which in linear response is given by: 

(s) = hTva (G A J X G R + G R J X G A ) AV E. (21) 

In s-space, each component is given as 

(s^ = hTTU^GiUGJ x GE. (22) 



We find (s z ) = (s x ) = 0, but 



The expression is simplified as 



-Jg^{6fc<? 1 + (A + A')<? }i?. (23) 

k 



s y ) = eAnrDXE, (24) 



by using the approximation (|19|). The spin accumulation is aligned to the pseudo-magnetic field of the spin-orbit 
interaction and its magnitude is proportional to the applied electric field within the linear response regime. The 
magnitude of the spin accumulation may be estimated as 



\(s y )\/D = 2.5 



A 



10- n eVm 



eE 



/' 



lfjScn^V^s- 1 



meV (25) 



10-KeV/m_ 

With A = 1CT 11 eVm0 the splitting of the chemical potential thus amounts to a significant 2.5 meV for typical 
experimental parameters of the applied field and the mobility /i. 

We have thus formulated the conductivity tensor and spin accumulation on an equal footing using linear response 
theory. The conductivity is found to be enhanced by the spin-orbit interaction and isotropic in spku space, i.e., 
all = ^xxi an d &Vx = = 0- We have proven that an electric field induced a spin-polarized densitycJ, but not a 
spin-polarized current. The recent discovery of the spin-galvanic effect, i.e. that a magnetization along y direction 
induces an electric currentO is reciprocal to the current induced spin accumulation. We may conclude that his current 
can also not be spin-polarized. 

The result that the conductivity is spin-isotropic implies that the spin accumulation in ferromagnet(F)/2DEG 
hybrids cannot be detected in two-terminal configuration with one ferromagnetic contact. A single source or drain 
ferromagnetic contact does not modify the global transport properties in the diffusive regime, because the contacts, 
which cpanect the reservoir distribution^fainctions to the semiconductor ones, are not affected by a magnetization 
reversalEZl. A phenomenological theoryEj'CJ is at odds with this conclusion. Microscopically, we trace the. matrix 
character of the current operatpij-as the culprit of this disagreementE3.__Experiments on_E/2DEG systemstil, which 
were supported by that theoryEfftj, were challenged by Monzon et aZo and van Weest3, who suspected that the 
measured effects were chie-to local Hall voltages caused by fringe fields near the ferromagnetic contacts. TheseEj and 
subsequent experimentsEj'Ej should perhaps be reconsidered in the light of the present theoretical results. 

The present results are related but different from the spin-Hall effect discussed by ZhangEl who dealt with a 
ferromagnetic metal thin-film. He suggested to measure the spin-accumulation excited by the electric current in a 
three-terminal configuration, however, which is also an option for the 2DEG, since we find the spin accumulation 
signal to be quite significant. 

Spin-polarized transport can be detected in a F/2DEG/F configuration with two ferromagnetic contacts as studied 
by Pareek and BrunoEa. We point out that the spin dependent conductances r^j and in Ref.Eil must not be 
confused with the spin dependent conductivities a x \ and aj, x defined above. The former conductances are defined 
as T^m) = (e 2 /fo)Trtf f (Dtj^-p where t is a transmittance matrix. Whereas &l x = as shown above, in general 
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rVi 7^ 0. It is possible to decompose the conductivity <t|J. by "cutting lines" in the conductivity diagram such that 
&xx = **xV^ + ^\x^ i w hich is simplified by the spin-diagonal vertex function in the ladder approximation. The two 
components are then given as 

hp 2 . 

&xx = Tc 7? + sb\k)G s (G + + G_) (26) 



and 

r TliT _ 



(7 



\QttL 2 ^ 

k,s 



hp 2 . 

J2 <b 2 k 2 + sbXk)G s (G+ - G_). (27) 



16ttL 2 

k, s 



We did not find a simple relation as Eq. (20) but in general cr]*P ^ o\NT , because in the limit A = 0, ffjj. = 0, 
whereas a^I tends to the Drude conductivity. This conclusion appears to be at odds with the numerical findings of 
Pareek and Bruno that in the limit of long samples ~ r^. 

In summary, we derived explicit expressions for the conductivity tensors and spin accumulation of a Rashba 2DEG 
with isotropic scattering centers, taking into account the vertex correction in the ladder approximation. The diffusive 
conductivity limited by non- magnetic impurity scattering is not spin dependent, although the applied bias does excite 
a spin accumulation. 
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